Universality of liquid- gas Mott transitions at finite temperatures 
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We explain in a consistent manner the set of seemingly conflicting experiments on the finite temperature 
Mott critical point, and demonstrate that the Mott transition is in the Ising universality class. We show that, 
even though the thermodynamic behavior of the system near such critical point is described by an Ising order 
parameter, the global conductivity can depend on other singular observables and, in particular, on the energy 
density. Finally, we show that in the presence of weak disorder the dimensionality of the system has crucial 
effects on the size of the critical region that is probed experimentally. 



Although band theory predicts that a system of electrons in 
a solid with one electron per site (unit cell) should be metal- 
lic, such a system ultimately insulates ma once the local 
electron repulsive interactions exceeds a critical value. The 
onset of the insulating state, the Mott transition, arises from 
the relative energy cost of the on-site Coulomb repulsion U 
between two electrons on the same lattice site, and the kinetic 
energy, represented by the band width W . Then, the transi- 
tion is governed solely by the ratio of U /W . At T = 0, it 
is often the case that symmetries of the microscopic system, 
associated with charge, orbital or spin order, may be broken in 
the Mott insulating state. However, at sufficiently high tem- 
peratures T, or in strongly frustrated systems, no symmetry is 
broken at the finite-T Mott transition. Then, the transition is 
characterized by paramagnetic insulating and metallic phases, 
whose coexistence terminates at a second-order critical point, 
depicted in Fig. |l(a)| In this paper, we are concerned with 
the universal properties of this classical critical point as 
revealed by a series of apparently conflicting experiments on 
(Cri_3;V;i;)203 101 and organic salts of the k - ET family 11]. 

Since no symmetry is broken at the finite-T Mott transi- 
tion, in a strict sense there is no order parameter. Nonetheless, 
experimental flSl. as well as theoretical evidence fllrtl sug- 
gest that the transition is in the Ising universality class, sim- 
ilar to the liquid-vapor transition. For example, Castellani et 
al. 101 constructed an effective Hamiltonian for this problem, 
and proposed that double occupancy should play the role of 
an order parameter for the Mott transition. On the insulating 
side, doubly occupied sites are effectively localized, but in 
the metal, they proliferate. A Landau-Ginzburg analysis 
provided further evidence for a non-analyticity in the double 
occupancy at a critical value ofU/W that defines a Mott tran- 
sition. Ising universality follows immediately because double 
occupancy, {fii-^riii), is a scalar local density field. 

Experimentally, the universality of the Mott critical point 
is typically probed by some external parameter, such as pres- 
sure, which can tune the ratio W/U . Measurements of the 
conductivity, S, on (Gri-xN x)20j found that away from 



the critical point, the exponents defined through 

AI](t,/i = 0) = = 0) - Sc oc It]''-, 

AS(t = 0,/i) oc \h\^'^\ 
d^{t,h)/dh\h^o oc (1) 

have mean-field Ising values, (3a — 1/2, 73- — 1 and 5a — 3. 
Here, t = {T - Tc)/Tc and h ^ {P - Pc)/Pc, with 
(Sc, Tc, Pc) denoting the corresponding values at the critical 
endpoint. Close to the critical region, Limelette et al. ^ ob- 
served a drift to the critical exponents of the 3D Ising univer- 
sality class. Mean field behavior is also seen in NiS2 iS]. 

However, similar pressure measurements ^ on the quasi- 
2D organic salts of the k-ET family appear to challenge the 
view that the Mott transition is in the Ising universality class. 
In this material, Kagawa et al. lH found that their data is de- 
scribed by the exponents (3^ — 1, 7(t — 1, and Sa — 2, 
which do not seem to be consistent with the known expo- 
nents of the 2D Ising model whose exponents are ||9|l /? = |, 
7 = 1 and (5=15. Since the exponents obey the scaling law 
ja = (3a {5a ~ 1), it was proposed that the Mott transition is 
in a new, as yet unknown universality class. The situation is 
further complicated by thermal expansion measurements ifioll 
that claim to measure the heat capacity exponent a and find 
0.8 < a < 0.95. This result is not only in sharp contrast to the 
expectation for an Ising transition (where q = for d = 2), 
it also strongly violates the scaling law a + 2(3a ^a = % if 
one uses the exponents of Ref. Jst]. 

In this paper, we present a unified phenomenological de- 
scription of all of these experimental facts within an Ising- 
type model, and resolve the issue of the universality class 
of the Mott transition. A complete description of these ex- 
periments requires to take into account that the conductivity 
depends on all possible singular observables of the associ- 
ated critical system, and not just on the thermodynamic or- 
der parameter associated with the phase transition. Similar 
considerations were made in magnetic systems near the Curie 
temperature ll ll Il2l 11311 . to explain the critical exponent of 
the conductivity along the coexistence curve. In that case, a 
symmetry of the microscopic definition of the conductance 
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FIG. I : (a) Typical phase diagram of Mott transitions as a function 
of pressure P and temperature T. At low T and P, a Neel antiferro- 
magnetic insulator (AFI) appears; it becomes a paramagnetic insula- 
tor (PI) if T increases, or a paramagnetic metal (PM) if P increases. 
Dashed line: first-order transition ending at a liquid-gas critical point. 
Full lines: continuous phase transition to the ordered state. Colored 
regions: critical (dark) and mean-field (light) regimes of the critical 
point, (b) An Ising configuration on the triangular lattice. Up (down) 
spins correspond to conducting (insulating) grains of linear size of 
the order of the system's dephasing length. 



prevented any coupling of the global conductivity to odd mo- 
ments of the order parameter, along the coexistence line. Even 
though similar in spirit, the situation here is much different. 
Starting from an effective microscopic model near an Ising 
critical point, we show that: 1) the conductivity typically de- 
pends on all possible singular thermodynamic observables of 
the system, namely the order parameter and energy density 
of the Ising model; 2) when the coupling to the energy den- 
sity dominates, there exists a large regime around the critical 
point, where the critical exponents for the conductivity are 
{(3a,Jcr,Sa-) ~ (1, |, ^), that agree (within the error bars) 
with the findings of Kagawa et al. ISJ, and the correspond- 
ing mean-field exponents are 7^^^, (5^^^) (1, i, f ); 
3) a crossover to Ising exponents is obtained in the order pa- 
rameter dominated regime as seen in Refs. |0, IM]; 4) in the 
presence of disorder the Mott critical point ultimately belongs 
to the random-field Ising model universality class, and there- 
fore the dimensionality of the system under study is even more 
important for specifying its critical properties. 

In order to resolve the discrepancies raised by these ex- 
periments, we consider the behavior of the conductivity of 
the system near the Mott critical point, assuming that it be- 
longs to the 2D Ising universality class. Rather than starting 
from a microscopic picture, e.g. a Hubbard model, we con- 
sider a coarse-grained model with the correct symmetries in 
which the physics of the relevant transport degrees of free- 
dom is captured. In this picture, one defines coarse-grained 
regions, of linear size of the order of the dephasing length 
^0 of the system, which are either insulating or conducting. 
Along these lines, we consider an Ising model on a 2D lattice 
(cf. Fig. |l(b)[ i. Near the critical point, where the correlation 
length for density fluctuations ^ diverges, it is expected that 
the relevant degrees of freedom behave classically. The Ising 
variables Si on each lattice site represent the fluctuating den- 
sity of mobile carriers on microscopic "grains" of linear size 
of the order of the dephasing length l^, which are conducting 



{Si = +1), or insulating {Si = —1). The Hamiltonian is 



where T is the temperature, P and Pc are the pressure and the 
critical pressure, respectively, and h oc P — Pc plays the role 
of the Ising magnetic field. This model is expected to describe 
the physics near the critical point, where ^ ^ Z^. In this limit, 
all other interactions beyond nearest-neighbor are irrelevant. 
Near the critical point, the most singular effect of the pressure 
is described by a coupling to the order parameter. 

To relate the order parameter fluctuations to the transport 
properties we will define an associated resistor network for 
this model, an approach that has been successfully used in 
other strongly correlated systems [Tsll . Let ac and ct/ be 
the local conductivities of the conducting and insulating re- 
gions, respectively. We define the bond conductance of the 
network model simply by adding these two conductivities in 
series. The bond conductance has three possible values, de- 
pending on the state of each grain, which can both be conduct- 
ing, both insulating, one conducting and the other insulating. 
Thus, the conductance of the bond has the form 



Cij = CTo (1 + gm{Si + Sj) + geStSj) 



(3) 



Even in this toy model, the microscopic conductivity, cTy , 
couples both to the order parameter. Si, and to the energy 
density, SiSj, of the Ising model with naturally large cou- 
plings, and g^, defined in Eq.Q. More specifically, we 
find that ao = i(ac + cti) + .g,„ = ^q^, and 

^'^'^""'^^ At high contrast, ac > 



9'^ 4(Jo{ac+(Ti) 

g,n — .9e — 1, whereas, at low contrast, \ai 



4(70 

<Ti, we get 



we get g^< gm^ 0. 

The conductivity of the 2D Ising model we described is 
a non-trivial quantity to compute. As it was shown in the 



simpler case of the random resistor network (RRN) Ill6[l . net- 
works of bonds with conductance ac (aj) chosen randomly 
with probability p and 1 —p, the global conductivity becomes 
non-zero as soon as an infinite percolating and conducting 
cluster emerges in the system. When aj = 0, the critical 
exponent (3a- of the conductivity is non-trivially related to the 
fractal properties of the incipient infinite conducting cluster 
This exponent is larger than unity for random uncorrected 
networks and larger than the exponent of the order parame- 
ter, because dangling bonds of the infinite cluster do not con- 
tribute to the conductivity. On the other hand, it becomes 
much less than unity for correlated networks, and typically 
very close to the exponent of the order parameter, since the in- 
finite cluster is efficiently connected with few dangling bonds. 

On the other hand, when ct/ > 0, a conducting cluster is 
less distinguishable from that of an insulating one, and the 
complex effects coming from the fractal cluster boundaries are 
smeared out. In the context of RRN, the percolation transition 
is not seen in the behavior of the conductivity, which seems 
to show just a crossover. If the contrast is low, aj c± ac, 
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FIG. 2: Crossover behavior of the conductivity at finite contrast: the 
energy density (order parameter) dominates at short (long) length 
scales. Inset: Fractal scaling at large contrast. 



FIG. 3: Monte Carlo data which verify the expected behavior of the 
conductivity when gm,ge ^ 1- E'^^°" scales as the energy density 
(see text). Inset: scales as the order parameter 



the actual conductivity of a single bond between sites E, 
should depend only on local observables, and we can formally 
expand it in powers of g,„ and lll7ll . 



O-Q + 9ni{{Si + Sj)) + fJeiSiSj) + 



(4) 



where the ellipsis represents more complex products of local 
spin operators (weighed by rapidly decaying functions) UtIi . 
Near the Ising critical point, the most singular contribution 
of the expectation values of multi-spin operators in Eq. |4] is 
given by the expectation value of the most singular, "primary", 
operators of the Ising critical point, the order parameter m and 
the energy density e. Thus, the most singular term of multi- 
spin operators with odd (even) number of spins is proportional 
to the order parameter (energy density). Therefore, within the 
range of convergence of this expansion, 

E ^ Eo(g™,5e) + fm{9m,ge){m) + f^{gm,9t){e) , (5) 

where 'Eq, fm, fe are non-universal regular polynomials in g,„ 
and g^. Provided that the critical behavior is still controlled by 
the fixed point theory of the Ising model, the total conductivity 
should have the structure of Eq. (|5]l. Thus, at finite contrast, 
Eq. Q predicts that the actual conductivity is the sum of even 
and odd components, under the Ising symmetry transforma- 
tion, E = Eo + E'=™" + E°'^^, and it should exhibit a crossover 
from an energy density dominated behavior at short distances 
to an order parameter dominated behavior at long distances. 
The crossover scale is controlled by the relative size of the 
functions /„i and (cf. Fig.|2]i. This behavior breaks down 
at high contrast where there is multi-fractal behavior (cf. Inset 
in Fig|2]and Ref. 111). 

We can understand the experiments of Refs. 10. IH S], if 
we assume that Eq. ^ applies. The results of Refs. 01 
follow by assuming that fmigm,ge) > fe{gm,ge), and the 
conductivity scales as the order parameter. Conversely, the 
results of Ref. 13] follow if fe{gm, ge) > fmig-m, .9e), and the 
conductivity, for an extended regime near the critical point, 
scales as the energy density of the Ising model. In this case 



holds AE oc |7Ti| , where 6 = (1 — a) Then, it follows 
that f]^ = op, 5a = 5/9 and 7,, = 7 + - 9). The 
resulting critical exponents are {[3a, ^a, 5a) = (1: |i "^)' very 
close to the experimental values. These exponents obey 



(6) 



if 7 = /? (5 — 1), i.e., the conductivity exponents obey a scal- 
ing relation identical to the Ising exponents, in agreement 
with the experimental verification of this scaling relation in 
Ref. ||5|l . In addition, the scaling function obtained by Kagawa 
et al. ISll only depends on [35 = [3a5a, as in our theory. 

In order to verify the theoretical picture presented above, 
we performed Monte Carlo simulations of the 2D Ising model 
on square and triangular lattices, using the Wolff cluster algo- 
rithm ifigll . For the calculation of the conductivity, for each 
Ising configuration we used the Franck - Lobb algorithm lEoll . 
or explicitly solved Kirchhoff equations. As expected, we 
found that at the Ising critical point, for gm,ge ^ 1, the even 
component of the conductivity E**™"^ scales as the energy den- 
sity, while the odd component E°'^'^ scales as the order param- 
eter (cf. Fig. [3]). As gm,gt approach unity, a slow crossover 
exists to a fractal regime of the Ising clusters, which is cru- 
cial for specifying the critical exponent of the conductivity, 
consistent with the results of Ref. llSll (cf. Inset in FiglH) 

Refs. 10, lit] report 3D mean-field Ising behavior and a small 
critical region in (Cri_2:Va;)203 and NiS2 underpressure re- 
spectively, in contrast to the extended critical region with 2D 
Ising exponents of Ref. We can understand these exper- 
iments by considering the effects of quenched disorder on an 
Ising critical point. The difference between a quasi-2D and a 
3D material is a strongly anisotropic Ising interaction along 
the direction perpendicular to the planes. Disorder that lo- 
cally favors the localized over the delocalized state or vice 
versa, corresponds to a random magnetic field, and couples to 
the order parameter of the Ising transition. This induces den- 
sity fluctuations. The relevant model for this discussion is the 
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anisotropic 3D random-field Ising model (RFIM), 



H 



E 



SiSj — Jz ^ SkSi + ^ 
{ki}. i 



h.Si 



(7) 



where hi is a random field with variance A. For = 3, there 
is a continuous phase transition in the 3D random-field Ising 
model (3DRFIM) universality class 112 ill for any anisotropy 
Jxy/Jz, an irrelevant operator at the 3D RFIM fixed point. 
However, for large anisotropy and weak disorder, relevant to 
the quasi-2D organics, there is a large dimensional crossover 
regime from 2D RFIM behavior, with an exponentially long 
correlation length, to the narrower 3D RFIM criticality 122 1. 
What changes between the 3D isotropic materials and the 
quasi-2D organics is not the universality class, but where 
the planar correlations become critical. For weak disorder 
A ^ Jxy and strong anisotropy Jz/Jxy ^ 1, the planes 
are essentially decoupled and 2D-RFIM behavior holds with 
£,xy ^ 1 in a large region away from the transition point. 
When Jz — Jxy, the critical region is narrow, and controlled 
by the 3D RFIM fixed point. 

With regards to the thermal expansion measurements that 
claim to measure the heat capacity exponent a, we argue that 
the authors of Ref. ifioll misinterpret their results. The volume 
change is proportional to the Ising order parameter of the Mott 
transition, i.e. Al cx m, yielding I'^dl/dT cx t'^^^. The 
thermal expansion diverges with exponent 1 — [3 = 0.875, 
consistent with Ref. ifioll who find it in the range 0.8 — 0.95. 

Some major predictions can be drawn from our picture. 
Firstly, all thermodynamic observables near the Mott criti- 
cal point should have Ising critical exponents. Secondly, re- 
garding the critical behavior in quasi-2D organic salts [51], in 
the clean system, the conductivity along the coexistence line 
should have the same critical exponent (/3cr = 1) in both mean- 
field and true-critical regimes. This means that the conductiv- 
ity jump AEj = S(T, h = 0+) - S(T, h = 0") along the 
coexistence line, which should be proportional to the order pa- 
rameter, should have distinct mean-field and critical regimes, 
where /^aSj = 1/2 and /JaSj =1/8 respectively. Also, 
the first-order Mott transition is expected to be broadened by 
disorder |23, 2^. Thus, instead of a sharp jump in the con- 
ductance one should see a continuous change which would 
become more abrupt for clean systems. The net effect is to 
make the system spatially inhomogeneous, as in charge or- 
dered phases, stripes and electron nematics 12511 and in the 
manganites l26ll . which tend to round their phase transitions 
and replace the first-order transition by an inhomogeneous 
phase. Thus, hysteretic glassy-like aging effects 12711 . a prob- 
lem that has been studied only recently in strongly correlated 
systems lfl4ll28ll . are also expected at this Mott transition. 

In conclusion, we have explained under a consistent phe- 
nomenological framework the series of experiments that were 
performed during the last few years on Mott criticality. We 
showed that the conductivity of a system near a critical point 
depends on all possible local singular observables of the sys- 
tem, which in the case of interest, are the order parameter and 



the energy density of the effective Ising model. This descrip- 
tion holds when the contrast between conducting and insulat- 
ing regions is small. Should this not hold, fractal behavior of 
the incipient infinite clusters at the critical point is crucial for 
specifying the critical exponents of the conductivity. Finally, 
disorder affects the effective dimensionality of the system. In 
particular, we showed that for quasi-2D materials, such as the 
organic salts in the k-ET family, critical fluctuations are ex- 
pected to be much larger than in a 3D material in an extended 
regime in the {P, T) plane around the Mott critical point. 
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